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We consider small network models for mutually delay-coupled systems which typically do not 
exhibit stable isochronally synchronized solutions. We show analytically and numerically that for 
certain coupling architectures which involve delayed self feedback to the nodes, the oscillators become 
isochronally synchronized. Applications are shown for both incoherent pump coupled lasers and 
spatio-temporal coupled fiber ring lasers. 
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I. INTRODUCTION 

Synchronization of networked, or coupled, systems has 
been examined for large networks of identical [1] and 
heterogeneous oscillators [2). For coupled systems with 
smaller numbers of oscillators, several new dynamical 
phenomena have been observed, including generalized [3[, 
phase [4], and lag Q synchronization. Lag synchroniza- 
tion, in which there is a phase shift between observed 
signals, is one of the routes to complete synchrony as cou- 
pling is increased [5[ and may occur without the presence 
of delay in the coupling terms. 

For systems with delayed coupling, a time lag between 
the oscillators is typically observed, with a leading time 
series followed by a lagging one. Such lagged systems are 
said to exhibit achronal synchronization. In [6], the ex- 
istence of achronal synchronization in a mutually delay- 
coupled semiconductor laser system was shown exper- 
imentally, and in [7], studied theoretically in a single- 
mode semiconductor laser model. In the case of short 
coupling delay for unidirectionally coupled systems, an- 
ticipatory synchronization occurs when a response in a 
system's state is not replicated simultaneously but in- 
stead is anticipated by the response system [H, Q, and 
an example of anticipation in synchronization is found 
in coupled semiconductor lasers [TO]. Cross-correlation 
statistics between the two intensities showed clear max- 
ima at delay times consisting of the difference between 
the feedback and the coupling delay. Anticipatory re- 
sponses in the presence of stochastic effects have been 
observed in models of excitable media as well as 
in experiments of coupled semi-conductor lasers in a 
transmitter-receiver configuration [29] . When the zero 
lag state is unstable and achronal synchroniza- 
tion occurs, the situation may be further compli- 
cated by switching between leader and follower. 
Switching has been observed theoretically and ex- 
perimentally in stochastic systems [12] but may 
occur even in deterministic chaotic systems [30]. 

Given that both lag and anticipatory dynamics may 
be observed in delay-coupled systems, it is natural to 
ask whether the isochronal, or zero lag, state, in which 
there is no phase difference in the synchronized time se- 
ries, may be stabilized in coupled systems. Stabilizing 
the isochronal state is important in bi-directional 
chaotic communication systems, as shown in re- 



cent theoretical work on communicating in sys- 
tems with delay [31]. Stable isochronal synchroniza- 
tion of semiconductor lasers has been observed recently 
in experiments [l3|, 03] and numerically [HI, EH- ^x- 
amples of partial isochronal synchrony in which 
only some of the oscillators in a delay network 
synchronize may be found in [27|, [3l|, and re- 
cently a theoretical explanation for partial syn- 
chronization has appeared in [28]. Other examples 
of isochronal synchrony have appeared in neural network 
models with delay [16|, [l7| . 

In this letter, we explore the possibility of adding self 
feedback to two globally coupled situations: 1. Inco- 
herent delay-coupled semiconductor systems [l8|, and 2. 
Coupled spatio-temporal systems consisting of coupled 
fiber ring lasers pj| with delay [20( | . 

We consider N coupled oscillators of the following 
form. Let F denote an m-dimensional vector field, B 
an m x m matrix, and Kj, where j = 1 • • • N, denote the 
coupling constants. For the cases we examine here, we 
consider global coupling including self feedback: 



dxjjt) 
dt 



F(Xi(t),Xi(t-T))+ E KjB Xj (t-T). (1) 



Given the structure of Eq. [TJ we examine the stability 
transverse to the synchronized state, S = {xi(t) : Xi(t) = 
s(i), i = 1 • • • N}, by defining rjij = X j —Xi. The linearized 
variations in the direction transverse to S are then given 

by 



*lhM=D 1 F(x i (t),x i (t-T))r lij (t) 
+D 2 F(xi(t),Xi(t - T))rHj(t - t) 
+(k» - Kj)Bxi(t — r) - KjBrjij(t - r) 



(2) 



where Di denotes the partial derivative with respect to 
the I th argument. 

We make the following hypotheses to simplify the anal- 
ysis: (HI): Assume that the dependence on the time de- 
layed variables in F takes the same form as the delay 
coupling; i.e., D2F(x,y) = Bnf. (H2): Let Ki = Kf = 
ft, i = 1 • • • N. Equation [2] then simplifies to 



drjjj (t) 
dt 



D 1 F(x l (t),x l (t-r))r ]ij (t), 



(3) 



where it is understood the arguments of the derivatives 
are computed along the synchronized solution s(t), and 
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the solution is a function of parameters such as coupling 
and delay. Computing Eq. [3] along the synchronized state 
will generate the Lyapunov exponents for the transverse 
directions, and we examine the effect of coupling and 
delay by computing the cross-correlations between time 
series as well. 

To examine the stability of the isochronally synchro- 
nized state of Eq.HJ we model N = 3 lasers that are pump 
coupled [H, [2l|. An isolated semiconductor laser's dy- 
namics at the i th node is governed by ^ = F(zi),Zi = 
(x h yi), where 



(a) 



F(z) = [— y — ex(a + by), x(l + ; 



(4) 



and x,y are the scaled carrier fluctuation number and 
normalized intensity fluctuations about steady state zero, 
respectively, e 2 is the ratio of photon to carrier lifetimes, 
and a and b are dimensionless constants (see [22| for de- 
tails on the derivation). 

The coupling strengths are K{ = ft/ = ft, z = 1,2,3. 
This leads to the following set of differential equations 
for the system: 



dzi(t) 
dt 



kJ2 B ^(t- 

i = 1 



r),i=l,2,3, (5) 



where m = 2, and 5(1,2) = 1, with all other entries 
in B equal to 0. An example of the intensities with 
and without self feedback in Fig. [T] shows explicitly the 
effect of self feedback in stabilizing the isochronal solu- 
tion. Writing down the differential equation for the trans- 
verse directions in matrix form for Eq. [5] using Eq. [3] 
and expanding near the synchronized solution rjij = 0, 
we obtain X'(t) = A(t, ft, r, e)X(t), where A(t, ft, r, e) = 
DF(s(t, ft, r, e)), and X(0) = I. Due to the nature of the 
global coupling with self feedback, each node receives the 
same signal. Therefore, the transverse stability does not 
explicitly depend on the coupling or delay, but rather on 
the dynamics of local nodes [23| • To examine the stability 
of the isochronal state, we derive some properties of the 
transverse Lyapunov exponents (TLE). The TLE satisfy 
the following limit: \(x ,y ,u) = lim t ->oo ~ log ^ffl" 11 ■ 
Here u is a vector in a given direction. 

By computing the solution to the linear variational 
equations along a given solution, we can extract the 
TLE. To examine the scaling behavior of the TLE, let 
A(t, ft, r, e) = det(X(t, ft, r, e)). Then, we have that 

A(t, ft, r, e) = exp( / trace(A(s, ft, r, e))ds) [24[. Taking 

the log of the matrix solution, and noting the determi- 
nant of a matrix is the product of its eigenvalues, we 
have: 

m -y 

VA(x ,2/o,ei) = lim - log | det(X(t, ft, r, e)|, (6) 

z — ' £^oo t 

i=l 

where are arbitrary independent basis vectors. Equa- 
tion [6] yields a rate of volume change in the dynamics 
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Figure 1: An example of delay-coupled dynamics showing in- 
tensities computed for N = 3, ft = 3.0c, t = 30, a = 2,6 = 
l,e = VO.OOl, using Eq. [4] (a) shows a solution where the 
lasers are coupled globally without self feedback, in which 
isochronal synchrony does not occur, (b) shows a stable 
isochronal solution with self feedback terms included. 



in the transverse directions. The solution may still be 
chaotic with one or more exponents being positive, but 
if sufficiently dissipative, volumes will shrink over time. 

From Eq. [H since trace(A(t, ft, r, e)) = — e(a + 
by(t, ft, t, e)) + x{t, ft, r, e), and assuming the inver- 
sion, x(t, ft, r, e), has zero time average due to sym- 
metry (which is observed numerically [26[), we have 
J trace(A(s, ft, r, e))ds = —e(a-\-b(y KyTy6 ))t and from 
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Figure 2: Prediction of the scaling of the sum of transverse 
Lyapunov exponents for Eq.[5]with respect to e. Other param- 
eter values are as in Figure [TJb) . Squares are the prediction 
using Eq. [T) and dots are the numerical values. 



Eq. [H we have 

AOo, 2/o, ei) + \(x , y , e 2 ) = -e(a + b (z/«, T , c ))- (7) 

Since e appears explicitly, it is easy to see how the sum 
of the TLE scales with e and compares with numerical 
experiments as in Fig. [2l 

Although the sum of the TLE is negative, loss of syn- 
chrony due to instability may occur at intermediate val- 
ues of e, as seen in Fig. [3j Regions where the isochronally 
synchronized solution is unstable are associated with one 
or more positive transverse Lyapunov exponenents. On 
the other hand, for sufficiently large damping, the trans- 
verse exponents reveal a stronger overall reduction in 
the phase space volume. The stability of isochronal syn- 
chrony with respect to other parameters can also be com- 
puted, e.g., as shown in Fig. [H for variations in coupling 
strength k. 

To illustrate the robustness of the self feedback struc- 
ture for generating isochronal synchronization in delay 
coupled systems, we examine a spatio-temporal stochas- 
tic system with multiple delays composed of coupled fiber 
ring lasers. A fiber ring laser system without self feed- 
back was studied in [20[ , and we extend the same model 
to include self feedback terms. In each ring laser, light 
circulates through a ring of optical fiber, at least part 
of which is doped for stimulated emission. The time for 
light to circulate through the ring is the cavity round-trip 
time tr = 202 ns, and the delay time in the coupling and 
self feedback lines is a second delay = 45 ns. Each laser 
is characterized by a total population inversion W(t) and 
an electric field E(t). The equations for the model dy- 
namics of the j th laser are as follows: 



Ej(t) = Re^[T(l-iaj)Wj(t)+i/±(l)}Ef h (t) 



dt 



q-l- Wj (t) 

-|4 db (t)| 2 {exp[2r^(t)]-l}. 



(8) 



(9) 




Figure 3: (Color online) (a) Transverse Lyapunov exponents 
and (b) cross-correlation (CC) of the dynamics for the same 
conditions as in Fig. [2] In (b), the cross-correlations be- 
tween lasers 1 and 2 (solid line) and 1 and 3 (dashed line) are 
shown. For most values of e shown here, a cross-correlation of 
1 is achieved when the shift between the time traces is zero, 
showing that the isochronal solution is stable. 




Figure 4: Cross-correlation (CC) between lasers 1 and 2 (solid 
line) and between 1 and 3 (dashed line) vs. coupling k for 
Eq. [5] Other parameters are the same as in Figure [TJb) . 
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The electric field from earlier times which affects the field 
at time t is 

Ef h (t) = Ej(t - Tr) + *l E l(t ~ Td) + KfEj(t - T d ). 

(10) 

Ej (t) is the complex envelope of the electric field in laser 
j, measured at a given reference point inside the cav- 
ity. E^ d (t) is a feedback term that includes optical feed- 
back within laser j and optical coupling with the other 
laser. Time is dimensionless. Energy input is given by 
the pump parameter q. Each electric field is perturbed 
by independent complex Gaussian noise sources, £j, with 
standard deviation D. We use a fixed input strength for 
all coupling terms: Ki = Kj = k for all i. (Values of the 
parameters in the model as well as further computational 
details can be found in [2(j. The only difference in pa- 
rameters was that the lasers are not detuned relative to 
each other in the current work.) 

Because of the feedback term Ej dh (t) in Eqs. El one 
can think of Eqs. [8] as mapping the electric field on the 
time interval [t — tr, t] to the time interval [t, t + tr] in 
the absence of coupling (k = 0). Equivalently, because 
the light is traveling around the cavity, Eqs. [8] maps the 
electric field at all points in the ring at time t to the elec- 
tric field at all points in the ring at time t + tr. We can 
thus construct spatio-temporal plots for E(t) or the in- 
tensity I{t) = \E(t)\ 2 by unwrapping E{t) into segments 
of length tr. 

Figure [5] shows time traces of the N = 2 lasers for 
a single round trip for both the system with self feed- 
back described here and the system without self feed- 
back (Kf =0) [25[. Isochronal synchrony can been seen 
when self feedback is included, while in the absence of self 
feedback the lasers are delay synchronized. The spatio- 
temporal plots in [5](c) and (d) are nearly identical due 
to the isochronal synchrony. To quantify the synchrony, 
we align the time traces for the two lasers with various 
time shifts between them. In the absence of self feed- 
back, the peak cross-correlation occurs when the lasers 
are shifted relative to each other by the delay time. The 
cross-correlation is low when the lasers are compared 
with no time shift. In contrast, when self feedback is 
included, the lasers achieve a high degree of correlation 
when compared isochronally. For the time traces shown 
in Fig. 03(a), the peak cross-correlation of 0.9554 occurs 
when there is no time shift, although the cross-correlation 
when shifted by the delay time is near as high (0.9549). 

We have swept the coupling strength n for the sys- 
tem of two lasers with self feedback and computed 
the average cross-correlation when the lasers are com- 
pared isochronally. Figure [6] shows that the lasers are 
well synchronized for input strengths as small as 0.1%. 
Isochronal synchronization can be produced when 
the lasers are detuned as in [20], but this requires 
stronger coupling and self feedback (not shown). 

For N = 3 fiber ring lasers, we have done a similar 
computation for cases with and without self feedback 
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Figure 5: (Color online) Intensity (arbitrary units) for two 
lasers coupled with k = 0.009. The left panels are intensity 
vs. time for laser 1 (bottom curve), and for laser 2 (top curve): 
(a) With self feedback, (b) without self feedback. Spatio- 
temporal plots corresponding to coupling with self feedback 
for (c) laser 1 and (d) laser 2. 
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Figure 6: Average cross-correlation vs. coupling for two cou- 
pled lasers with self feedback. 



(not shown). We found that when the the lasers are cou- 



pled globally without self feedback, the isochronal state 
will still synchronize. However, adding self feedback will 
cause the isochronal state to stabilize at somewhat lower 
values of coupling. Further details for this case are in 




In summary, we have considered delay-coupled systems 
and, through the addition of self feedback, obtained sta- 
ble isochronal synchrony in coupled semiconductor and 
fiber ring laser models. Model analysis for incoherent 
pump coupled lasers reveals scaling of the Lyapunov ex- 
ponents transverse to the synchronized state, while com- 
putations on systems of coupled fiber ring lasers show 
how self feedback may cause the onset of synchrony in 
coupled spatio-temporal systems. In the cases we have 
studied, we constructed small globally coupled networks. 
For the small clusters presented here with delay, it is ad- 
vantageous to add feedback loops, since this was key to 
stabilizing the synchronous state. A question for future 
study is how this method may be scaled up for larger 
networks. 
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